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ABSTRACT 

Autofocus ttclmlquts or* Ming dttlgntd at tht «)tt Propulsion 
. Laboratory to autoMtIcally choose tha filter parameters (1.e.» the 
focus) for the digital synthetic aperture radar correlator; currently, 
processing relies upon Interaction with a hunan operator mho uses his 
subjective assessment of the quality of the processed SAR data. 
Algorithms have been devised applying luge cross-correlation to 
aid In the choice of filter parameters, but this method also has Its 
drawbacks In that the cross-correlation result may not be readily 
Interpretable. Enhanced performance of the cross-correlation tech- 
niques of JPL was hypothesized given that the Images to be cross- 
correlated were first filtered to Improve the signal-to-noise ratio 
for the pair of scenes. The results of experiments are described 
and images are shown. 
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1.0 INTRODUCTION 


Tht Jtt Propulsion Libontory ntlnulns • digital synthttic aptr- 
turt radar corrtlator «N1ch Is now bting utlllztd to process the raw 
data fron the Seasat«A SAR alsslon which took place In 1978. The data 
represent the coinplex video signal* and they swst be filtered to pro- 
duce Images suitable for viewing. Currently* filter parameters are 
chosen by a human operator* based solely upon subjective evaluation of 
Image ’’quality". Alternative means for determining the "focus" para- 
meters haw been sought by jPL. One technique consists of cross- 
correlating lamges generated as, say* the first and fourth looks which 
originate from distinct spectral components of the complex SAR Images. 

As discussed In JPL Internal memoranda [1* 2]* when the SAR data 
are not optimally focused, the first and fourth look Images will be 
offset from each other spatially, and the offset can, theoretically, 
be determined through registration techniques (using cross-correlation 
algorithms). However* the Images generated have statistically Inde- 
pendent noise degradation (I.e.* independent fading noise) since the 
scenes arise from different portions of the complex spectrum of the 
SAR data. Although the underlying microwave reflectivity scene Is 
Ideally the same for the Image being correlated, a strong correla- 
tion peak cannot usually result due to the formation of different 
noise sample functions In the two scenes. 

It was proposed that the cross-correlation approach would work 
better If the two scenes were first filtered to Improve signal-to-nolse 
characteristics, prior to crors-correlatlng. An adaptive minimum mean 
square error Image filter developed by Frost et al. [3] at the Remote 
Sensing Laboratory was candidate for preprocessing two one-look Images 


3 


"amplitude" Image ims squared so that Intensities could be displayed. 
Next, the Image was segmented Into "homogeneous" regions as determined 
by field boundaries (the scene Is an agriculturally developed region 
* in Iowa). 

The noisy areas within the field boimdarles were "colored over" 
on the digital Image display at the Remote Sensing Laboratory, with 
the mean values of those areas determining the Intensity that would 
replace the fading scene. Next, this segmented, averaged Image (shown 
In Figure 2) was Input to an algorithm which regenerates fading values 
according to the relation [5. 6] 

<’> 

where N represents the number of looks for the radar (assumed Integer) 

and y Is a random variable whose probability density function Is given 

2 

as a X 2 f)‘ Thus, a random number generator was employed to replace 
average Intensities. No pixel correlation was used to color the spec* 
trum of the noise; this approach Is supported by the observation that 
the pixel-to-pixel correlation Is not highly significant for N ■ 1 
Images even though there are several pixels per resolution distance 
(actual experiments were conducted with the SAR 580 N • 1 data, which 
was digitally processed, and not sampled at the large Interval (* 17 m) 
of the commercially available Seasat*A SAR data). Since white noise 
was used. It Is possible the results of the overall experiment would 
change somewhat If samples were correlated (say, run through a linear 
filter prior to multiplication with the "defaded" scene). 

Two N « 1 Images were computer generated with different noise 
sample functions as seen In Figure 2.. These Images were next Input 
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The cross'correlatlw surface peak and Its surrounding region for 
Images 1 and 2 Is shonn In Figure 6. It Is characterized by a weak 
correlation maxiaun, and It Is pyramidal In shape, with no strong side- 
lobe structure. This surface was photographed from an Isometric dis- 
play; a 64 by 64 pixel region from the original 512 by 512 pixel sur- 
face Is represented. The equivalent rectangular width of this surface 
Is 53 pixels by 53 pixels. The peak value of the cross-correlation Is 
listed In Table 4. The cross-correlation of Images 3 and 4 Is viewed 
In Figure 7. The adaptive filtering effects have caused the cross- 
correlation peak to have more sharply crested zones away from the peak. 
The equivalent rectangular width Is 57 pixels square; the cross-cor- 
relation peak Is listed In Table 4. Likewise, the cross-correlation 
of Images 5 and 6 (processed with Identical filter parameters (B) of 
Table 4) Is shown In Figure 8. Again, It Is noted that a more sharply 
crested fali-off of the peak Is produced by the WISE filtering. The 
equivalent rectangular width Is 59 pixels square; the cross-correlation 
peak Is listed In Table 4. 

The three cross-correlations are shown together In Figures 9 and 
10. Note that all Images were scaled Identically so that a direct 
comparison could be made. To examine more closely the fine structure 
of the peaks, they were rephotographed as side views, and they are 
shown In Figures 11, 12, and 13. We see that although the cross- 
correlation of Images 1 and 2 (original fading Images), which Is 
given as Figure 11, has a smaller equivalent rectangular width than 
the preprocessed cross- correlations. It Is apparent that the surface 
is sharper near the peak for Figures 12 and 13 (Image 3 cross-corre- 
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estimator and research Is continuing along these lines. 

Previous research of Fries and Nodes tino [8] has shown the statis- 
tical limitation of the autocorrelation peak for scenes similar to the 
"source” Image of Figure 2. Thus, It is known that the cross-correla- 
tion peak cannot rise above a certain level on the average, even If 
the two Image Inputs to the cross-correlation consisted of Identical 
values r(t) on a pixel -by-pixel basis due to preprocessing. That Is, 
the type of scene limits Its own unnormallzed autocorrelation peak. 
Fries and Nodestino also have shown this effect to be functionally 
dependent on edge density and average Intensity correlation between 
neighboring fields [8]. One could consider the agricultural test 
region chosen to be somewhat a worst case example of SAR scenes for 
analysis. The greatest difference In cross-correlation peaks for 
nonprocessed versus preprocessed scenes would be seen for two Inde- 
pendent, white noise fields; this case might be thought of as the best 
case for demonstrating potential of the MMSE filtering approach. 
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Figure 1. Multiplicative noise model for radar signals followed 
by cross-correlation, and then by minimum mean square 
error filtering as a preprocessing step for cross- 
correlation. 
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Output 


Figure 3. Cross-correlation surface generation. 
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Figure 5. 


Image 5 was generated from Image 1 using the adaptive 
filtering parameters "B" of Table 3. Image 6 was com 
puted from Image 2, also based upon parameters “B" of 
Table 3. 
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Figure 7. 




Pictured above is the cross-correlation surface for 
Images 3 and 4. Cresting of the surface is seen; the 
shape is still primarily pyramidal. 
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Figure 12. Side view 
of cross-correlation sur- 
face of Images 3 and 4. 


Figure 13. Side view 


of cross-correlation sur 
face of Images 5 and 6. 
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'^^Surell. Side view 
- of cross-correlation sur- 


face of Images 1 and 2. 
Identical scales were 
used for photographs on 
this oage. 
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Table 2 

Filter Paraneters and Thei-" Functions 


Parameters 


Functions 


Filter Size 

Number of Filters 
(NF) 

Observation Window 
Size 

Maximum Alpha 

Input Gain 

Alpha Bias 
Alpha Gain 


filter impulse response wjdth 

sets NF-1 thresholds for which 

are then mapped into filter indices 

2 2 

window size in which o and u are 
computed 

sets the rate of decay within a filter 
as a function of filter index 

multiplying factor on the input pixel 
values 

adjusts the filter indices up or down 

maps the observed noise-to-signal 
ratio (o2/p2) »nto a filter index 
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Table 4 

Cross-Correlation Results 


Image Number 

Correlation Peak 
U 10^ 

Equivalent 
Rectangular Width 

1 and 4 

1.648 

53.0 

2 and 5 

2.096 

58.7 

3 and 6 

2.020 

58.4 

7 (autocorrelation) 

1.986 
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The minimum mean square error filter for stochastic estimation 
of fading radar signals based upon the multiplicative SAR image model 
[5, 6] was developed and reported previously [7, 3]. The implemen- 
tation of this filter, in effect, determines whether a region under 
observation is "homogeneous" (e.g., a white noise field) or not (con- 
tains a boundary, or for some other reason possesses a signal-to-noise 

ratio lower than the theoretically calculable value). NF-1 thresholds 

2 2 

exist with which to compare the observed value of o /w , where 


(ws)" 

(ws) 


(ws)‘ 


T*; 5'"- 


ws equals the window dimension (square in this implementation). 

2 2 

The quantized levels of a /u are mapped into filter indices, 

NF total filters. The filter chosen is used as a mask to multiply on 
a pixel -by-pixel basis with the underlying image data. As the obser- 
vation window moves, different filters are picked to apply to the 
data. Within a filter, the exponential decay of weightings is governed 
by the constant alpha max. 

The following computer output lists the filters used on Image 3 
(generated from Image 1) and gives the number of filter usages. 

The adaptive implementation of the MMSE filter w,s written by 

* 

Stephan A. Smith and Victor Frost of the Remote Sensing Laboratory. 


* 

Now at Texas Instruments, Dallas, Texas. 
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